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HIGH- CURRENT  BETATRON,  PHASE  I 


1-1  INTRODUCTION 

In  a  conventional  3etatron,  the  current  is  limited  by 
the  method  of  injection  and  the  constant  orbit  radius.  For 
example,  with  an  injection  energy  of  100  keV  and  an  orbit 
radius  of  1  m,  the  initial  magnetic  field  must  be  ^10  G.  The 
current  density  is  limited  by  space  charge  of  the  beam  to 
'vlOO  mA/cm^.  We  proposed  (MLR-402-1,  April  28,  1980)  to  con¬ 
struct  a  device  with  a  toroidal  magnetic  field  of  ^10^  G,  in 
which  case  the  space  charge  limit  for  the  current  density  is 
increased  by  the  factor  10^.  We  proposed  to  inject  the  elec¬ 
trons  to  form  a  beam  by  an  adaptation  of  the  method  of  induc¬ 
tive  changing. 

1-2  STATEMENT  OF  WORK  (MLP-402-1) 

The  Statement  of  Work  from  proposal  MLP-402-1  is 
reproduced  below. 

1-2.1  Statement  of  Work  (November  1980-June  15,  1981) 

Task  1.  Construct  Torus 

A  glass  torus  with  a  40  cm  major  radius  and  5  cm 
minor  radius  will  be  constructed,  coated  with  epoxy 
fiberglass  and  wound  with  copper  wire. 

Task  2.  Design  and  Fabricate  Coils 
Vertical-field,  air-core  coils  will  be  designed  and 
fabricated. 

Task  3.  Assemble  Systems 

The  systems  of  the  Betatron  experiment,  including 
magnetic  field  coils,  power  supplies,  vacuum  system, 
and  injectors,  will  be  assembled. 


Task  4.  Perform  Initial  Experiments 
Initial  experiments  on  injection  trapping  and  acceler¬ 
ation  will  be  performed.  Electron  current  and  energy 
will  be  measured. 

Task  5.  Analyze  Stability 

Electron  orbit  and  beam  stability  will  be  investigated 
analytically  and  numerically.  Stability  analyses 
will  be  performed  for: 

a.  Integer  and  half  integer  resonances. 

b.  Negative  mass  instability. 

c.  Resistive  wall/negative  energy  wave  insta¬ 
bilities  . 

This  work  was  carried  out  during  the  period  November 
1980- June  15,  1981.  The  results  are  described  in  Appendix  A. 

1-2.2  Statement  of  Work  (June  15-December  31,  1981) 

Phase  I  was  continued  during  the  period  June  15- 
December  31,  1981.  The  Statement  of  Work  for  this  period  is 
reproduced  below: 

Task  1.  Measurement  of  Injected  Charge  and 
Trapping  Efficiency 

The  injected  charge  can  be  determined  from  current 
traces.  The  charge  trapped  can  be  measured  by  observ¬ 
ing  the  frequency  of  diocotron  oscillations.  The 
charge  will  first  be  trapped  in  local  mirror  regions 
as  illustrated  in  Figure  1.  Then,  as  the  diamagnetic 
current  in  the  copper  washer  decays,  electrons  will 
be  ejected  along  the  toroidal  field  lines.  It  is 
important  to  measure  the  electron  line  density  before, 
during  and  after  this  process.  Parameter  studies 
will  be  carried  out  including  variations  in  injector 
voltage  and  injector  position.  We  plan  to  begin  with 
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Figure  1.  Betatron  experiment. 


one  injector  and  increase  the  number  up  to  four,  as 
required.  The  end  result  of  this  study  should  be 
that  a  toroidal  beam  of  electrons  is  trapped  and  has 
stable  confinement  for  'vl  ms.  The  mean  velocity  in 
the  toroidal  direction  should  be  zero. 

Task  2.  Acceleration  of  the  Toroidal  Electron  Beam 
The  Betatron  fields  will  be  applied.  This  involves  a 
peak  magnetic  field  of  1  kG  at  the  location  of  the 
beam  with  a  rise  time  of  ^1  ms.  The  toroidal  electric 
field  will  produce  ^100  V  per  revolution.  The  field 
index  can  be  varied  from  0.2  to  0.8.  We  will  measure 
the  current  with  a  Rogowski  coil  and  monitor  the 
charge  per  unit  length  by  continuously  observing  the 
diocotron  frequency.  X-ray  measurements  with  scintil¬ 
lators  will  also  be  used  to  monitor  the  beam.  The 
objective  is  to  document  acceleration  to  a  few  MeV 
with  a  current  of  the  order  of  1  kA. 

1-3  SUMMARY  OF  RESULTS  FOR  THE  PERIOD 

JUNE  15 -DECEMBER  15,  1981 

A  local  mirror  was  produced  around  the  electron 
injectors  by  weakening  the  toroidal  magnetic  field  with  a 
closed  loop  of  wires.  Electrons  were  trapped  there  as  well  as 
in  a  driven  mirror.  The  line  density  of  trapped  electrons  was 
proportional  to  the  injection  voltage.  With  the  use  of  four 
injectors,  electrons  were  captured  in  four  local  mirrors  along 
the  torus.  The  current  in  the  loop  of  wires  decayed  much 
faster  than  the  toroidal  field  resulting  in  collapse  of  the 
local  mirror.  Electrons  were  then  lost  due  to  toroidal  drift 
as  indicated  by  a  delayed  burst  of  X-rays. 

1-4  EXPERIMENTAL  PROGRESS 

Figure  2  shows  schematically  how  the  apparatus  has 
been  modified.  In  a  previous  report,  an  arrangement  was  des¬ 
cribed  whereby  the  mirror  field  was  produced  by  energizing 


mirror  coils  with  a  capacitor  bank  (Figure  2a) .  This  approach 
was  subsequently  abandoned  and  a  passive  means  was  adopted.  A 
pair  of  stranded  cables,  2.5  mm  dia,  was  wound  around  the 
injector  (Figure  2b).  These  cables  were  shorted  to  form  a 
closed  loop.  An  electric  current  was  induced  by  the  shorted 
cables  by  the  changing  toroidal  magnetic  field,  producing  a 
local  mirror. 

After  electrons  were  successfully  trapped  in  the 
mirror  formed  by  the  shorted  cables,  two  additional  injectors 
were  installed  and  tested  (Figure  2c).  The  multiple  injectors 
were  first  fired  independently,  the  A-K  gap  was  adjusted,  and 
some  conditions  (e.g.,  heater  current)  for  electron  trapping 
were  checked.  Electrostatic  probes  were  located  in  only  one 
mirror.  For  the  remaining  two  mirrors,  the  anodes  of  the 
electron  injectors  were  used  as  elec'  ostatic  probes.  After 
each  injector  was  tested  independv  all  three  were  fired 

simultaneously  and  electrons  were  trapr  id  in  three  mirrors. 

The  experimental  apparatus  was  disassembled  to  install  additional 
electrostatic  probes  and  an  additional  injector.  The  final 
configuration  involved  four  injectors,  seven  shorted  cable 
pairs  and  seven  electrostatic  probes  (Figure  2d) . 

1-5  INDUCTIVE  CHARGING  IN  A  PASSIVE  MIRROR 

A  comparison  between  the  driven  mirror  and  the  passive 
mirror  is  made  in  Figure  3.  In  a  driven  mirror,  the  electro¬ 
static  oscillation  lasted  longer  (Figure  3a) .  In  the  passive 
mirror,  the  oscillation  damped  in  150  ys  (Figure  3b).  The 
oscilloscope  traces  of  the  oscillations  were,  however,  quite 
similar  to  each  other  when  observed  on  a  fast  sweep  speed. 

This  suggests  that  the  trapping  takes  place  similarly  in  both 
cases . 

Figure  4  shows  the  electrostatic  oscillations  detected 
by  anodes  when  three  injectors  were  used.  Simultaneous  trapping 
was  observed  as  in  Figure  5.  The  frequencies  of  the  oscillations 
in  the  three  cells  were  not  identical  but  were  nearly  the  same. 


vV 


Figure  4. 


MLI82-409 


Oscilloscope  traces  of  electrostatic  oscillations 
picked  up  by  the  anode  of  electron  injector. 
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Figure  5.  Electrostatic  signals  from  two  mirrors.  50  ys 
Upper:  6  «  0°.  Lower:  9  *  135°. 


The  dependence  of  trapping  on  experimental  parameters 
was  studied.  Figure  6  shows  oscilloscope  traces  of  electrostatic 
oscillations  when  changing  the  injection  voltage.  All  the 
other  parameters  were  fixed.  The  frequency  is  proportional  to 
the  injection  voltage  (Figure  7)  in  agreement  with  theory. 

The  frequency  is  plotted  as  a  function  of  injection 
current  in  Figure  8.  It  is  noted  that  the  frequency  is  weakly 
dependent  on  the  injection  current:  the  frequency  variation  is 
within  10  percent  for  the  total  injection  current  of  0.4-2. 0  A. 
The  trapping  was  sensitive  to  the  injection  timing,  as  studied 
on  the  collective  focusing  ion  accelerator. 

1-6  LIFETIME  OF  ELECTRONS 

The  mirror  ratio  changes  with  time.  Figure  9a  shows 
the  current  which  flows  through  one  section  of  the  toroidal 
coils.  Figure  9b  is  the  current  in  one  of  the  closed  loops. 

These  waveforms  indicate,  roughly  speaking,  that  the  mirror 
ratio  decreases  rapidly  after  the  time  when  the  toroidal  field 
reached  its  peak.  An  X-ray  burst  appears  soon  after  this  time 
as  shown  in  Figure  9c.  Electrostatic  oscillations  were  not 
observed  afterwards. 

1-7  FOUR  MIRROR  EXPERIMENT 

The  injectors  used  in  the  four  mirror  experiment  were 
not  all  of  the  same  design.  Two  of  the  injectors  (0  *  0°, 

135°)  gave  electron  trapping  in  wide  ranges  of  the  experimental 
parameters,  but  the  other  two  (0  -  180°,  315°)  worked  in  very 
limited  ranges.  This  made  the  simultaneous  trapping  in  four 
mirrors  difficult  and  unrepeatable.  Insufficient  time  and 
resources  were  available  to  fabricate  new  injectors  of  a  common 
design. 

In  addition  to  electrostatic  probe  measurements  in 
the  mirror  cells,  the  electrostatic  oscillations  were  detected 
by  an  electrostatic  probe  located  halfway  between  0  *  135°  and 
180°.  The  frequency  measured  between  the  cells  was  the  same  as 
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Figure  9.  Waveforms  from  toroidal  experiment. 
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that  measured  at  0  *  180  within  the  limits  of  experimental 
error  and  repeatability. 

1-8  DISCUSSION 

The  process  of  injection  has  been  verified  for  a 
small  number  of  mirrors.  The  present  experiment  is  limited  to 
magnetic  fields  between  3. 5-7.0  kG  and  a  small  number  of 
inj ect or /mirror s .  Electrons  are  trapped  and  contained  as 
expected  in  the  mirror  fields.  However,  when  the  mirror  fields 
are  relaxed  there  is  insufficient  charge  to  avoid  the  toroidal 
drift.  (This  is  consistent  with  estimates  of  the  line  density 
of  the  injected  charge  as  determined  from  the  diocotron  oscilla 
tions.)  As  a  result,  the  electrons  drifted  to  the  walls  when 
the  mirror  fields  were  relaxed  instead  of  being  contained  in 
the  torus. 

This  research  will  continue  under  ONR  sponsorship. 
Some  of  the  limitations  of  the  experimental  apparatus  relating 
to  power  conditioning,  coil  design,  etc.,  will  be  improved. 

Thus  far,  no  new  physics  problems  have  arisen  and  the  princi¬ 
ples  of  the  injection/ trapping  scheme  have  been  documented. 

A  continuous  theoretical  effort  has  been  carried  out. 
The  main  results  to  date  are  a  comprehensive  analysis  of  the 
problems  of  orbital  stability  that  are  presented  in  Appendix  B. 
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HIGH  CURRENT  BETATRON* 


G.  Barak,  0.  Chernon,  A.  Fisher,  H.  Ishizuka,  and  N.  Roscoker 
University  of  California,  Irvine,  California,  and 
Maxwell  Laboratories,  Inc.,  San  Diego,  California 

ABSTRACT 

In  a  conventional  Betatron,  space  charge  limits  the  electron  current  during 
injection.  By  adding  a  toroidal  magnetic  field  and  injecting  electrons 
with  the  inductive  charging  method  we  expect  to  increase  the  current  limit 
by  a  factor  of  10^  -  10".  After  the  beam  has  been  accelerated  to  high 
energy  the  toroidal  field  is  no  longer  necessary  and  can  be  reduced  to  fac¬ 
ilitate  extraction  of  the  beam.  The  effect  of  the  toroidal  magnetic  field 
on  space  charge  instabilities  is  investigated.  An  experimental  Betatron 
has  been  constructed  with  a  major  radius  of  40  cm  and  a  minor  radius  of  3 
cm.  On  the  basis  of  previous  experimental  results  with  inductive  charging 
the  current  limit  should  be  a  few  kilo-amperes. 

I.  INTRODUCTION 


In  a  conventional  Betatron  the  orbit  radius  is  R  ■  Ymc^/eB  where  B  is 

2  — k  7  7. 

the  Betatron  field,  3  *  v/c  and  Y  •  (1-9  ).  For  the  usual  injector1 

Y  •  1,  3  <  1;  for  example,  with  a  100  keV  injector  and  R  ■  1  meter,  B  » 

2  ^ 

10.6  gauss.  The  space  charge  current  limit  is  proportional  to  B^  ,  and  is 
therefore  very  small. 

The  current  limit  has  been  increased  by  reducing  the  orbit  radius  and  in- 

2 

creasing  the  injection  energy. .  Another  proposal  to  eliminate  the  space 
charge  problem  by  Introducing  a  plasma  was  not  successful.3  Controlling 
space  charge  of  electrons  by  adding  a  toroidal  magnetic  field  to  the  usual 

4 

Betatron  fields  was  first  investigated  theoretically  in  1973,  and  the 
first  experiments  were  reported3  in  1976.  A  similar  proposal  for  ions  was 
made  in  1978. 6 

During  the  past  five  years  the  technology  of  inductive  charging  has  been 

developed  at  U.C.I.  in  experiments  with  magnetic  mirrors7  and  a  bumpy  torus 

A  line  density  of  4  xlO11  electrons/cm  has  been  injected  and  confined  for 

8 

several  milliseconds.  A  new  Betatron  experiment  has  been  constructed  dur¬ 
ing  the  past  year.  It  is  illustrated  schematically  in  Fig.  1.  A  local 
mirror  is  created  into  which  injection  and  trapping  take  place,  exactly  as 
in  the  previously  studied  bumpy  torus.  However,  the  diamagnetic  current  in 
the  copper  ring  (washer)  decays  on  a  short  time  scale  compared  to  the  toroi 
dal  magnetic  field.  The  electrons  trapped  in  the  mirror  are  then  ejected 
along  the  field  lines  as  the  mirror  collapses.  3ased  on  the  results  with 
the  bumpy  torus  we  expect  to  trap  enough  electrons  to  produce  a  currant  of 


a  few  kilo-amperes  when  the  electrons  are  accelerated  to  9  ■  1. 

Various  instabilities  may  arise  in  the  modified  Betatron  including  space 

charge,  negative  mass, resistive  wall  and  orbital  resonance  instabilities. 

In  general  the  toroidal  magnetic  field  is  stabilizing.  In  this  paper  we 

shall  consider  only  space  charge  instabilities;  the_ other  instabilities 

9  10 

have  been  considered  in  previous  publications.  ’  Space  charge  instabili¬ 
ties  have  previously  been  discussed  for  the  field  index  s  -  %  which  is  a 
special  case  in  which  the  analysis  is  simple.  The  present  analysis  inclu¬ 
des  s  and  a  consideration  of  electric  and  magnetic  images. 

II.  SPACE  CHARGE  INSTABILITIES 


To  describe  the  orbits  in  a  toroidal  electron  beam  consider  the  local  co¬ 
ordinates  (x,y)  or  (r,3)  illustrated  in  Pig.  2.  The  equations  of  motion 

sx  sB_ 

for  standard  Betatron  fields  B  *  B  (1  +  r— ) ,  B  »  — —  v  and  a  toroidal 

q  7  0  K  X  K 

field  B  ,  are  as  follows9 
z 


*  Cl  B 

x+^i  +  nzy  +  C(l-s)ny2-  n2]  X  -  -  g^y  (1) 

z 

y  +  7  y  -  nzi-  [atiy2 -n2]  y  -  ^  ^x  (2) 

Z 


CJ  •  eB ,/Ymc,  ■  eB  /Ymc,  ^/2Y*  where  «  *  «  4-TT  ne^/Ym  and  s  is 

z  z  7  o  p  p 

the  field  index.  (x,y)  describe  oscillations  about  the  equilibrium  orbit 
given  by  the  Betatron  condition  R  ■  v^/Q^  where  s-c  is  the  velocity  of 
the  electron,  in  the  toroidal  direction. 

For  the  special  case  s  •  .5  these  equations  may  be  combined  via  the  trans¬ 
formation  ^  • 

S  -  (x+iy)  Y^exp[-^^  dt'r*2(t')]  (3) 

o 


to  give 


* 

* 


U) 


2. 


u> 


-I 


where 


-n2  + 


1  d 
1  Ht 


(4) 


Transverse  beam  stability  requires  tu  >  o.  For  a  coasting  beam  (Y 
corresponds  to 


v  2  12 

n  <®y  *,  > 


4rrmc 


o)  this 
(5) 


The  advantage  of  using  a  strong  toroidal  field  is  apparent  when  Y  **1  and  3y 
is  small.  After  acceleration  when  Y  and  By  are  large  the  toroidal  field 
can  be  permitted  to  decay  to  facilitate  beam  extraction. 

The  limiting  current  is  determined  by  the  injection  conditions.  For  a  con- 


(6) 


ventional  3etatron  a  <  YB^  /4rrmc  implies  a  maximum  current 


I  -  (17/4)  y3  3  3  (r.  /R)2  k-amps 

ZQ  D 


For  zz*  nodlfiad  Betatron  with  a  toroidal  field 


Im  -  (17/8)‘y3  33(rb/R)2(B2/B0)2  k-amps 


(7) 


where  we  have  used  the  Betatron  condition  R  *  3c/ft 


If  s  is  not  exactly  .5  the  general  solution  of  Eqs. (1)  and  (2)  gives  three 


independent  conditions  that  must  be  satisfied. 
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^2+fiz2) 


.22 

4<JU  'M 
x  y 


2  2 

“*  “y  > 


ui 


uj  2  +  n  2  > 
y  z 


(8) 

(9) 

(10) 


where  w 


(l-s)ft2-ft2,  «J2-  sft2  -  a2  and  w2  - 
7  7  7 


2.  2 
u)  -Hu 
X  V  "z 

2  4 


r. 


We  have 


assumed  s,  Y»  0  .  0  and  0  are  time  Independent.  Equations  (8)  and  (10) 

z  F 

are  satisfied  if  m2>o.  However  Eq.  (9)  is  independent  of  ft  .  This  con¬ 


dition  is  illustrated  in  Fig.  3  where  Q ft/fi^)  is  plotted  against  s  and  the 

stable  region  is  indicated.  In  the  lover  region  conventional  Betatrons 
2  2 

operate  and  w  ,  uu  >  o.  In  the  high  current  Betatron,  the  initial  operat¬ 


ing  point  is  in  the  upper  region,  but  as  Y  rises  we  eventually  must  cross 

the  unstable  region  except  for  s  ■  1/2.  In  the  unstable  region  the  beam 

will  expand,  ft2  will  drop  and  a  transition  to  the  lower  stable  region 

will  occur.  For  example,  if  s  ■  4(1+6)  the  density  drop  will  be  In/n  as 

2  I  Al 


If  N  is  the  line  density  and  r.  is  the  beam  radius  the  transition 
3  2  °  2 

occurs  for  Y  such  that  Y  3  >  4Nr  (R/r,  )  where  r  is  the  classical  elec- 

e  o  e 


tron  radius.  Assuming  N  ■  10^-  cm"^,  R  ■  40  cm,  r,  ■  1  a  the  transition 

D 


takes  place  at  Y  *  6 . 

The  Fermi  drift  has  not  been  included  in  Eqs.  (1)  and  (2).  The  Fermi  drift 


is  due  to  tht  fact  that  the  toroidal  field  is  not  constant;  i.e.,  B  *  B 

[l  +  so  that  the  drift  velocity  iA  y»v^/2,ftzR.  In  Eq.  (1)  ftzy  must 

be  replaced  by  ft.Cl  +  r]  y.  The  Fermi-drift  term  C5  x  y/R  is  non-linear. 
zo  *»  2  2 

It  may  be  compared  to  the  term  ft  x  or  (l-s)ft  x.  In  the  first  case  if 
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*»o  7/R 


ft- 


y  V2  F i\2 

"  “  S  (uT/  <<  1 


P 


the  Fermi  drift  term  can  be  neglected.  Physically  the  drift  is  corrected 
by  the  self-field  drift  which  causes  angular  rotation  at  the  frequency 


uj  «  n  /H  .  The  displacement  due  to  Fermi  drift  y/w  is  negligible  for  a 
o  zo  o 


dense  beam  and  small  Y  .  If  Y  is  large  the  self  electric  and  magnetic 
fields  almost  cancel.  However,  in  that  case  the  Fermi  drift  may  still  be 
neglected  if 


2f<  y/R 
zo  } 


a  2 


(  px  /  Ymc )  <<  1 


In  this  case  the  Fermi  drift  is  compensated  by  the  focusing  effect  of  the 
Betatron  field.  The  Fermi  drift  may  be  neglected  except  when  (l-s)ft 


uix  —0.  This  is  an  unstable  region  in  any  case  and  should  be  passed  quick¬ 


ly  as  discussed  above. 

If  the  electron  ring  is  in  a  conducting  torus,  but  displaced  from  the  axis 
electric  and  magnetic  images  at  the  wall  will  create  forces  on  the  beam 
that  cause  it  to  drift.  Equations  of  motion  for  the  center  of  mass  of 
the  beam  moving  in  external  fields  including  image  fields  are  similar  to 
Eqs.  (1)  and  (2).  The  results  are  formally  the  same  with  the  revised  de¬ 
finitions 

2 


(1-s)  n  2-  (rb/a)2  n2 


«y  -  %2  -  (rb/a)2n2 


n2  -  ujp2/2 


We  have  assumed  that  image  currents  in  the  wall  decay  immediately  and  image 
charges  survive,  r^  is  the  beam  radius  and  a  is  the  minor  radius  of  the 
torus. 

The  beam  will  drift  around  the  toroidal  axis.  As  Y  increases  it  must  pass 

2  2  11-1 

through  an  unstable  region  when  Y@  ■  4Nre(R/a)  .  For  N  ■  10  cm  ,  R  « 

40  cm  and  a  »  5  cm,  Y~8.  In  this  case  the  self  field  term  r.  ^n2/a2  de- 

D 

pends  only  on  line  density  which  is  not  altered  by  the  Instability.  It  is 
necessary  to  have  s^l/2  and  to  change  y  rapidly  enough  that  the  beam  does 
not  displace  to  the  wall. 

The  unstable  region  for  beam  motion  is  not  the  same  as  the  unstable  region 
for  single  particle  orbits.  In  the  neighborhood  of  the  latter  unstable 
region  the  Fermi  drift  is  significant.  The  beam  drift  due  to  electrostatic 
Image  forces  should  then  be  considered.  The  drift  frequency  is  <*>o  « 

(r  2y2/a2)  Cl2/D  where  C?  ■  \u  * tly" ,  The  Fermi  drift  may  be  neglected  if 

D  ZO  p 


III.  EXPERIMENTAL  PARAMETERS 


An  experiment  has  been  constructed  with  the  objective  of  producing  a  I  k- 

aapere  beam  of  electrons  at  10  MeV.  The  torus  is  a  glass  chamber  with  a 

metal  screen  liner  of  major  radius  40  cm  and  minor  radius  5  cm.  A  vacuum 
—8 

of  5 x  10  torr  has  been  achieved  with  a  cryo-pump .  The  toroidal  magnetic 
field  has  a  peak  value  of  12  k-Gauss  with  a  rise  time  of  100  usee.  After 
crow-barring  the  decay  time  is  about  300  usee.  The  Betatron  field  has  a 
peak  value  of  1  k-Gauss  at  the  majof  radius  with  a  rise  time  of  350  usee. 
The  field  index  varies  from  .2  to  .8  in  the  chamber.  Magnetic  field 
errors  have  been  determined  to  be  less  than  12.  Injectors  are  similar  to 

g 

those  previously  employed  in  the  bumpy  torus.  They  operate  at  15  k-volts 
with  a  current  at  5  amps,  for  5  usee.  They  have  been  operated  at  voltages 
as  high  as  30  k-volts.  Construction  and  field  measurements  were  completed 
in  May  1981  and  injection  experiments  are  now  in  progress.  Figure  4  is  a 
photograph  of  the  U.C.I.  Betatron  experiment. 

*Work  supported  by  the  Defense  Advance  Research  Project  Agency  and  the 
Office  of  Naval  Research. 

We  are  indebted  to  Donald  W.  Kerst  for  many  illuminating  discussions. 
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Figure  3.  Stability  condition  for  the  modified  Betatron. 

Initial  and  final  operating  points  are  indicated. 
The  .transition  takes  place  when  Y  is  increased 
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Abstract 

We  analyze  orbital  stability  of  the  High  Current  Betatron. 

In  this  modified  Betatron  a  toroidal  magnetic  field  is  added  to 
the  conventional  Betatron  magnetic  field.  This  increases  sub¬ 
stantially  the  space  charge  limit  during  injection.  It  gives  rise, 
however,  to  new  problems  during  acceleration,  such  as  ring  sta¬ 
bility,  Fermi  drift  and  orbital  resonances.  A  detailed  analysis 
is  presented  showing  that  two  serious  problems  may  arise; (i)  The 
beam  becomes  highly  unstable  when  the  net  focusing  force  on  it 
becomes  zero.  The  dominant  instability  in  this  region  is  the  Fermi 
drift.  This  may,  however,  be  utilized  in  extracting  the  beam. 

(li)  Errors  in  the  vertical  (Betatron)  magnetic  field  result  in 
driven  resonances  that  must  be  crossed  during  acceleration.  The 
amplitude  growth  of  the  transverse  Betatron  oscillations  due  to 
these  resoances  is  negligible  up  to  very  high  y  (low  n  resonances). 
It  may  become  practically  very  difficult  to  cross  the  n-1  reson¬ 
ance.  Thus  again  extracting  the  beam  by  making  use  of  the  Fermi 
drift  while  n>l  is  preferable. 
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I.  INTRODUCTION 

The  largest  conventional  Betatron  was  completed  about  1950.* 
The  energy  was  about  300  MeV  and  the  beam  current  about  100  milli- 
amperes .  The  maximum  beam  current  was  limited  by  space  charge  dur¬ 
ing  injection.  Electrons  from  a  100  KeV  electron  gun  were  injected 

into  an  orbit  of  radius  R  *•  1  meter  and  a  Betatron  magnetic  field 

o 

of  By  ••  10  Gauss.  The  space  charge  limit  is  proportional  to  By 

1  7  -3 

and  for  the  Kerst  Betatron  n  ^  10  cm  .  The  space  charge  limit 

is  also  proportional  to  y,  thus  it  is  greatly  increased  after 

acceleration,  in  a  conventional  Betatron  it  is  the  space  charge 

limit  at  injection  which  determines  the  maximum  current. 

To  eliminate  the  space  charge  problem  the  plasma  Betatron  was 
2  3 

proposed  by  Budker.  After  many  studies  were  carried  out,  the 
maximum  current  reached  about  10  amperes,  much  less  than  expected. 
The  current  is  probably  instability  limited.  The  precise  insta¬ 
bility  has  not  been  identified,  although  the  negative  mass  insta- 

4 

bility  is  mentioned  frequently. 

By  increasing  the  injection  energy  and  decreasing  the  orbit 

radius  in  a  conventional  Betatron,  the  space  charge  limit  can  be 
2  3  2 

Increased  (n«£  y  /R  ) .  Small  "ironless"  Betatrons  have  been  de¬ 
veloped  with  an  electron  energy  of  100  MeV  and  electron  currents 

5 

of  about  90  amperes.  Electrons  were  injected  into  a  3.9  cm  or¬ 
bit  radius  with  energy  of  500  KeV.  The  space  charge  limit  was 

4 

increased  by  a  factor  of  ~  10  compared  to  the  Kerst  Betatron. 

The  total  number  of  electrons  was  Increased  however,  only  by  a 
factor  of  ~  40. 


Controlling  space  charge  of  electrons  by  adding  a  toroidal 

magnetic  field  to  the  usual  Betatron  field  was  first  investigated 

theoretically**  in  1973  and  the  first  experiments  were  reported*0 

12 

in  1976.  A  similar  proposal  for  ions  was  made  in  1978. 

The  combined  toroidal  and  vertical  (Betatron)  fields  increase 

substantially  the  space  charge  limit  and  the  number  of  electrons 

that  can  be  trapped  and  accelerated.  Injecting  and  trapping  into 

such  combined  fields  require  new  techniques.  Several  methods 

6  7  6 

were  proposed  lately.  *  The  first  method  has  already  been 

demonstrated  in  the  studies  of  the  Collective  Focusing  Ion  Accel- 
8 

erator.  The  method  is  called  inductive  charging  and  was  first 
a 

used  in  HIPAC.  Electrons  are  injected  from  thermionic  injectors 
and  trapped  by  means  of  a  rising  toroidal  magnetic  field.  A 
Betatron  with  this  type  of  injection  is  illustrated  in  Figure  1. 
Electrons  move  in  an  orbit  around  the  torus.  They  are  not  drifted 
in  the  toroidal  field  due  to  the  self  field  forces  and  the  inter¬ 
action  with  the  surrounding  metallic  boundaries.  There  are  no 

restrictions  on  Bz  as  there  are  on  By  in  a  conventional  Betatron. 

*  2 

As  the  space  charge  limit  depends  on  B„  during  injection,  a  sub- 

z 

stantial  increase  in  the  density  can  be  achieved.  With  the  exist¬ 
ing  injectors,  the  current  after  acceleration  will  be  ~  1  KAmpere. 
Better  injectors  may  increase  this  current  further  by  1-2  orders 
of  magnitude. 

The  toroidal  magnetic  field  solves  the  problem  of  increasing 
the  space  charge  limit  during  injection.  It  creates,  however, 
new  problems  during  acceleration.  It  is  the  purpose  of  this  paper 
to  investigate  the  regions  of  instability  of  this  high  current 


Betatron  and  to  suggest  an  optimum  mode  of  operation.  We  shall 
assume  throughout  this  paper  that  inductive  charging  is  used. 

The  parameter  set  we  shall  use  as  a  standard  reference  is  given 
in  Table  1.  We  shall  discuss  only  orbital  instabilities.  Collec 
tive  instabilities,  such  as  negative  mass,  or  resistive  wall  in¬ 
stabilities  have  not  been  treated  in  this  paper,  but  it  is  known* 
that  toroidal  field  decreases  their  importance. 


II.  SINGLE  PARTICLE  ORBITS 


2.1  Equations  of  Motion 

In  order  to  describe  the  particle  orbit  in  a  toroidal  elec¬ 
tron  beam,  consider  the  local  coordinates  (x,y)  or  (r,6)  as  illus¬ 
trated  in  Fig.  2.  The  toroidal  direction  is  indicated  by  z.  If 

eBv  vz 

n  “  electrons  will  move  on  a  radius  R  -  •*—  .  Furthermore, 
y  Ymc  Uy 

if  By  is  changed  in  time  so  that 


By(R) 


1  CR 

-“v  \  B  (r) r  dr 
2R  °0  y 


the  electrons  will  move  on  a  constant  radius.  Due  to  the  beam 
self  fields  and  the  distribution  of  the  initial  conditions  elec¬ 
trons  will  also  move  in  the  transverse  directions.  This  motion 
determines  the  maximum  current  that  can  be  obtained  and  the  quality 
of  the  beam.  The  equations  of  motion,  up  to  first  order  in  (x/R) , 
where  x  is  the  deviation  from  the  Betatron  equilibrium  orbit,  are: 


+  *=x]  -  '  e[«r-8Be)f  +|  (VyB2-»zBy)] 

-y[y+iy]  -  -  etCEr-eB^f+I^B^-v^)] 


(2.1) 


Bz  and  By  are  the  applied  toroidal  and  vertical  fields,  Ep  and  B0 
are  the  self  fields.  Interaction  with  the  surrounding  wall  was 
neglected  here  and  will  be  considered  separately  later.  The  equili 
brium  orbit  is  determined  by  the  relation  By 00-  h  ByOO  •  ByOO  is 
influenced  by  the  beam  current.  For  a  high  current  Betatron  this 
is  not  negligible,  especially  for  the  early  stage  of  acceleration. 


External  currents  should  be  added  to  eliminate  this  effect.  We 


shall  assume  that  this  is  done  so  that  v  *  R  Cl  ,  where  R  is  con- 

z  y 

stant  and  vz  “*  c  (the  case  of  Y  **  1  will  be  investigated  separately). 

For  a  cylindrical  symmetric  beam  the  self-fields  are  easily 
calculated, 

-2  ,  , 

-e(Er-!Be)/.,r-^  • 

_ 2  4Trne2 

where  u)p  -  — ,  n  is  the  average  density  of  the  beam  and  n(r) 
is  the  average  density  up  to  radius  r.  If  the  beam  is  uniform, 
n(r)/n  -  1,  and  the  self  force  is  harmonic. 

The  Betatron  field  is 


B  -  B 
y  o 


(_s_y 

\R-x/ 


Bo  (1  +  ir) 


B. 


sBoy 


(2.2) 


Thus  we  get  the  following  equations  of  motion: 

x  +  x  i  +  nz  y  +  [(l-s)fiy2  -  n2]  x  -  0 

y  +  y  ^  -  fiz  x  +  [s  ny2  -  n2]  y  -  o 


eBz  ®Bo 

where  -  - —  ,  fl  -  —2 
z  vac  ’  y  ymc  ’ 


n 


. 

2Y 


(2.3) 


The  toroidal  field  was  assumed  constant.  Higher  order  corrections 
will  be  discussed  later. 


2.2  Space  Charge  Limits 


Define 


-  (1-s)  G 


2 

y 


G- 


(2.4; 


and  assume  that  sf  y,  Gz,  fiy>  fi  are  time  independent.  The  general 
solution  of  Eqs.  (2.3)  and  (2.4)  can  be  expressed  as  a  sum  of 
two  oscillations  having  frequencies 


U) 


1,2 


n  2  +  u a.  2  +  uu  2  ± 
z  x  y 


:°z2+%2+*v2>  -4^2. 


“2) 


1/2 


(2.s: 


Orbits  are  stable  if  and  only  if: 

(a)  (uu2  4  ou2  4  fi2)  >  4w  2m  2 

*  y  *  y 

2  2 

(b)  %  <«y  >  0  (2.6 

(c)  «>x2  +  ««y2  +  nz2  >  0 

Denoting  s  *  1/2(14-6)  these  conditions  can  be  written  as: 

(a)  (nz2  +  20y2-4D*)  nS!2  +  «2ny4>  o. 

2 

(b)  C(n/ny)2  - 1/2]  >  |  62  (2.7; 

(c)  nz2  ♦  ny2  -  2n2  >  o 

Condition  (b)  Involves  only  0  and  Gy.  For  6-0  it  is  always 
satisfied.  When  6*0  it  is  satisfied  in  two  disconnected  regions: 


-  -  -  I M 


(JL 

ny; 


1+  6 


(2.8) 


The  first  is  the  "low  density"  region  in  which  conventional  Beta- 

2  2 

trons  are  operated  («J  ,  a)  >  0) .  The  second  is  the  "high 

x  y 

2  2 

density"  region  (uu  ,  id  <  0).  In  the  high  current  Betatron  we 

x  y 

intend  to  start  in  the  high  density  region,  but  as  y  rises  we 

eventually  have  to  cross  the  unstable  region  (see  Fig.  3).  Here, 

2 

the  beam  will  start  to  expand,  but  as  it  expands  ft  will  drop 
and  a  transition  to  the  lower  stable  region  will  occur.  An  upper 
limit  estimate  of  the  density  drop  is  |<^H|  <  2  |  6  | ,  thus  for  |  6  |  « 
0.25  the  beam  will  expand  by  less  than  25%  in  radius. 

This  effect  can  be  decreased  by  crossing  the  unstable  region 
quick  enough;  e.g.,  by  a  rapid  change  of  y. 

Conditions  (a)  and  (c)  are  satisfied  for  every  6  if  they 
are  satisfied  for  6-0  (s-1/2).  Condition  (a)  is  then  stronger 


than  (c).  Thus  we  get: 


2  nz2 

ft  <  +  -J-  or 


n  <  y  [ +  -isj] 

4Ttmc  8nmc 


(2.9) 


For  the  inductive  charging  By  «  Bz  at  injection, 

2 

thus  n  „  «  y  B_  .  In  experiments  at  the  University  of  California, 

max  z 

Irvine,  electron  ring  of  line  density  ~  2  •  1011  cm”1  was  injected 
and  trapped.  This  corresponds  to  a  current  after  accelera- 

4 

tion  of  about  1  Kamp,  and  is  larger  by  a  factor  of  10 


compared  to  the  Kerst  Betatron.  This  is  still  two  orders  of  magni¬ 
tude  less  than  the  current  that  can  be  contained  by  a  10  kG  toroidal 
field.  We  expect  to  increase  the  density  of  the  trapped  electron 

ring  by  using  higher  injection  voltage  which  is  presently  10  K-volts 

12 

In  experiments  at  Maxwell  Laboratories  a  line  density  of  2.10  was 
observed  with  an  injector  voltage  of  30  K-volts.*0 

During  acceleration  YBy2  rises,  thus  if  Eq.  (2.9)  is  satisfied 
during  injection  it  will  be  also  satisfied  during  acceleration. 

The  transition  from  the  high  fi/Qy  region  to  the  low  fi/Gy 
region  occurs  at 


3  2  4  N  R  r 

Y  9  “  - nr- — ' 


(2.10) 


where  r  is  the  electron  classical  radius  and  r.  is  the  beam 
c  b 

radius.  For  our  standard  example  (Table  1)  this  means  Y  “*  13 


and  B 


210  Gauss. 


III.  RING  ORBITS 


Assume  that  the  electron  ring  is  shifted  away  from  the  axis  of 

the  torus.  Electric  and  magnetic  images  at  the  wall  will  push  the 

beam  further  onto  the  wall.  Together  with  the  toroidal  magnetic 

field  a  slow  rotation  (F  x  B  drift)  around  the  toroidal  axis  will 

start.  We  shall  investigate  here  this  transverse  motion  and  find 

its  stability  regions.  We  shall  consider  only  the  z-independent 

mode,  and  expand  the  image  forces  to  the  lowest  order  in  a/RQ 

(where  a  and  Rq  are  the  minor  and  major  radii  respectively). 

To  the  lowest  order  in  a/RQ,  the  image  forces  may  be  considered 

as  created  by  an  infinite  cylinder  of  radius  a.  If  the  beam  is  off 

center  by  6x  along  the  x  axis,  the  image  will  be  a  line  with  the 

2 

same  line  densityN  located  at  x«a  /6x.  The  electric  field  at  6x 
resulted  by  the  image  charge  is 


-2Ne6x 


aZ/6x  -  6x| 


(3.1) 


and  the  force  exerted  on  an  electron  at  6x  will  be: 


F 


Y  m  fJZ(rb/a)Z 


(3.2) 


where 

2  _  "p  _  2nne2 

n  ~2 - y® 

The  magnetic  field  depends  on  the  image  currents,  if  the  wall 

is  an  ideal  conductor  the  magnetic  force  is  opposite  to  the  elec- 

o 

trie  force  and  reduces  it  by  a  factor  of  y  •  We  shall  assume  the 


other  extreme  case;  i.e.,  immediate  decay.  The  equations  of  motion 
for  the  beam  center  are: 


11. 


"  Vi  ,  z 
x  +  Y  +  V* 


,/r.  ,2  eB 

^  '  °zy  +  vif  *= 


(3.3) 


;  .  2/rb\  .  eBx 

y  +  }y  -  n  (— )  y  +  nzx  -  vs 


x  and  y  are  measured  relative  to  the  toroidal  axis.  If  the  beam 

center  is  shifted  relative  to  this  axis,  the  equilibrium  radius 

will  be  determined  by  a  different  condition  than  that  of  a  single 

particle,  namely  R  *  v  /nw(R).  This  is  due  to  the  additional  force 

z  y 

in  the  x  direction.  Assume  that  the  equilibrium 

orbit,  R,  is  shifted  in  the  x  direction  by  6x;  i.e.,  R  »  Rq  -  fix. 

If  the  Betatron  field  is  expanded  around  R  we  get: 


1  +  s 


x-6x 


B  ¥ 
o  R 


Note  that  v^  -  Rfi^R)  and  s  -  s(R),  BQ  -  BQ(R). 

Substituting  in  Eq.  (3.3)  we  get  for  the  x  component 


(3.4) 


*  +  7*  *  nzy  +  [<i-ony2  -  Q)  n2]  X 


-  (s-l)  r,  6x  (3.5) 


If  we  make  the  transformation 


x  -  x  - 


(s-i)n„  fix 


(l-s)  Qy  -  (rb/a)V 


(3.6) 


which  is  permitted  as  long  as  the  denominator  is  not  zero,  the 
equations  of  motion  become: 


x+^x  +  Q  y  +  uj  2  x  “  0 

I  Z  X 


(3.7 


•*  Y  •  •  2 

y+^y  -  n2  x  +wy  y  -  o 


where 


9  9  2  9 

V  =  (i-s)ny  -  (Va)  n 


(3.8 


2  2  /  v2  2 
V  "  S  V'  <rb/a)  fi 


The  beam  equilibrium  radius  is  thus  shifted  by 


Ax 


(l-s)fi  6x 

- - k-7t  relative  to  the  toroidal  axis.  The 


(l-s)fiy  -  (rb/a)~fi 


physical  reason  for  this  is  that  the  image  force  supplies  part  of 
the  centripetal  force. 


Equations  (3.7)  are  similar  to  Eq.  (2.3).  For  s  -  1/2  the 
motion  is  stable  if 


*  <  (£) 


w  •  i: 


L  4 


or 


(3.9 


2  B  2 


2 


,  V  “I 

Y  Vr.  /  L“ - 2  +  --y  oJ 

b  8rrmc  4rrmc^ 


This  equation  can  be  written  approximately  as: 


13. 


4.10 


+  2.10s7  y 


where  N  is  the  line  density  and  the  parameters  of  Table  1 

12 

were  used.  If  N< 1.8.10  ,  the  beam  is  stable  for  every  y. 

12 

However,  for  N?  1.8.10  there  is  a  region  of  instability. 

13  -1 

For  N  » 10  cm  this  region  starts  already  at  y“"40.  To 

overcome  it  either  B  should  be  raised  or  the  currents  in 

z 

the  wall  should  persist  longer. 

Equations  (3.7)  have  the  same  stability  problem  as  we  saw  in 
Eq.  (2.8).  .  It  will  start  at 


02  *  \  (*/rb)2  ny2 


or  Yp2  «  4N(R/a)2r 


(3.10) 


where  r  is  the  electron  classical  radius.  For  our  standard  param- 
c 

eters  this  will  occur  at  y  «  90.  In  contrast  to  the  6ingle 

particle  case,  the  self  field  depends  here  only  on  the  line  density, 

which  does  not  change  as  the  instability  develops.  It  is  thus 

necessary  to  cross  this  region  fast  enough.  An  estimate  of  the 

2  1  2  2 

growth  rate  can  be  given  at  the  maximum;  i.e.,  at  fi  (a/rb)  . 
Here  the  e-folding  time  is 


2nz  M  0.045 

1—41 


lisec 


which  is  much  less  than  nsec  unless  |s-$|  «  1. 

Besides  the  exponential  growing  instability  the  beam  will 
drift  in  the  y  direction.  This  is  due  to  the  r.h.s.  term  of  Eq. 

0  2 

(3.5).  The  drift  speed  is  j  6x  «  20  6x  cm/iasec.  Assuming 

6x  -  0.5  cm  the  beam  will  hit  the  wall  in  0.5  (isec.  Increasing 
2 

NR  will  increase  the  Y  where  the  transition  occurs.  Thus  the 


instability  may  be  avoided  (if  NR  is  large  enough)  or  it  may  be 

used  as  a  means  to  extract  the  beam. 

In  concluding  this  section  a  remark  should  be  made  concerning 

the  single  particle  orbits.  In  Eqs .  (2.1)  we  neglected  the  force 

created  by  the  image  charges  This  force  will  add  oscillating 

terms  to  Eqs.  (2.3).  A  problem  will  arise  only  when  the  single 

electron  and  the  beam  motions  are  in  resonance;  i.e.,  uu±(single 

particle)  -  uu±(beam).  This  is  not  likely  to  happen  because  the 
2  2 

u>  ,  a)  in  Eqs.  (2.4)  and  (3.8)  are  different  (as  long  as  fi  is 
a  y 

not  very  small  compared  to  fi  ). 


IV.  BEAM  EMITTANCE 


The  quality  of  a  beam  is  sometimes  described  by  its  emit- 
tance,  defined  as: 

€/tt  -  e  p-pA  P/pz  (4.1) 

where  p  is  the  beam  radius  and  8  is  the  angular  divergence.  Emit- 
tance  is  determined  mainly  by  the  method  of  injection.  For  in¬ 
ductive  charging  p  x  p  —  m  v^o  a  -  const,  where  a  5  cm  (the  beam 
radius  during  injection),  and  vi0/c«  0.2  corresponding  to  10  keV 
Injector  .  When  the  toroidal  field  reaches  its  maximum  of  ~  10 

KGauss  the  beam  radius  is  of  the  order  of  1  cm  and  the  transverse 

8 

kinetic  energy  is  ~500keV.'  During  acceleration  p  -  const,  but 
p  decreases  due  to  the  By  increase  and  the  sel f-pinching  effect. 

For  simplicity  assume  s  -  £  and  consider  Eq.  (2.3).  If  B_  is  chang 

ing  in  time,  the  induced  field  Eg  should  be  included.  Trans¬ 
forming  then  to  the  rotating  Larmor  frame  we  get: 


u  +  u i  +  £  (  nz2  +  2  Qy2  -  4  fi2  )u  «  o 
where  (4.2) 

x  +  iy  -  expr  £  (t')  dt'  j  u  (t) 

This  is  correct  for  any  change  of  fi  :  either  changing  Y  or  B. 

Z  Z 

.  v 

The  mixed  term  u  can  be  eliminated  by  the  transformation 

u  -  Y~^  z  (4.3) 


yielding 

where 


..  2 
Z  +  U)  2  «■  O 


+  2  Oy2  -  4  n2 


*unz2  +  2ny2-  4f2]  (4‘4) 


ii 


Eq.  (4.3),  for  the  real  and  imaginary  parts  of  z  can  be 
solved  using  the  WKB  approximation.  The  general  solution  can 


16. 


be  written  as 


ci  /  \ 

zi  -  j=- cos  r* +  v 

where  c^  and  cpi  are  determined  by  the  initial  conditions,  thus 
r2  -  y”1  U!2-  [a2  +  B2  sin(2u)t  +  cp) J 

2  2 

where  A  ,  B  and  cp  are  functions  of  the  initial  conditions  only. 
The  maximum  value  of  r,  thus  the  beam  radius,  changes  in  time  as 
(yu))“^f  or: 


r(t)  _  u(t)  .  v~^(t)  z(t) 
rtoT  utoT  z(o5 


(nz2  +  2ay2-4n2)o  (4.5) 

At)  2^-4°^' 


During  inductive  charging  “■  0  and  y  “•  1.  If  ft  “Owe  get  the 
well-known  adiabadic  condition 


r  (t)  B  (t )  -  const 

V  2 

This  is  equivalent  to  ■  const. 


(4.6) 


As  acceleration  takes  place  Qy  rises  and  Q  drops  down 

o 

(like  y  ) .  At  high  y  the  self  term  is  always  negligible  due  to 
cancellation  of  the  electric  and  magnetic  forces.  The  beam  rad¬ 


ius  then  changes  as 


r(t)  (Bz2  +  2  By2)  ^  -  const 


(4.7) 
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Thus  the  beam  is  further  compressed  as  it  is  accelerated. 


If, after  B  reaches  its  maximum  3Z  is  turned  off,  the  final  beam 
y  i 

radius  will  be  rf/rQ  “  (B^/2By2f)  .  For  -  10  kG,  Byf  -  5  kG, 


rf/ro 


1.2. 


At  a  final  energy  of  300  MeV  €/tt 


v,  1.2r 
xo  o 


Yc 


2  millirad  cm. 


Instabilities  may  increase  this  number.  As  the  emittance  de¬ 
pends  mainly  on  the  injection  method,  a  full  understanding  of 
this  process  will  be  necessary  to  further  increase  the  beam 
quality. 

The  electron  ring  dynamics  during  acceleration  is  very  simi 
lar  to  the  one  electron  dynamics.  Here,  however,  the  electro¬ 
static  interaction  with  the  walls  does  not  reduce  during  acceler 
ation  and  according  to  (4.5)  the  amplitude  Increases  as: 


r(t) 


2  B_ 

1  +  - fc-  -  2  y 

V 


m 


* 


-  const 


(4.8) 


For  the  standard  parameters  the  change  in  r  is  negligible, 
even  for  y  “  600.  If,  however,  the  line  density  is  Increased  by 
an  order  of  magnitude,  the  beam  will  become  unstable  at  y  «  200. 
This  can  be  remedied  by  increasing  the  toroidal  field  or  by  the 
persisting  image  currents  in  the  walls  (compare  to  Eq.  (3.9)). 


V.  FERMI  DRIFT 


A  charged  particle  moving  alo;_  .oroidal  magnetic  field 
drifts  with  a  constant  speed  in  a  direction  perpendicular  to  the 
toroidal  plane.  This  drift  is  known  as  the  Fermi  (sometimes  also 
VB)  dirft.  It  is  due  to  the  non-constancy  of  the  magnetic  field 
(necessary  for  v x B  -  0)  as  a  function  of  the  radius.  This 
effect  is  of  second  order  in  (a/R)  and  was  neglected  in  the  previous 
discussion.  In  the  Betatron  we  describe,  which  has  both  toroidal 
and  vertical  magnetic  fields,  and  which  has  strong  self-fields, 
the  Fermi  drift  motion  is  weakened  as  we  shall  see. 

Let  us  start  with  the  simple  case  of  electron  moving  along 
a  toroidal  fields.  The  equations  of  motion  in  the  transverse 
directions  are: 


where 


x  +  Q„(x) y  -  0 

y  -  nz(*>  i  -  o 

-  R 

”  (r-x)  nz(1+  R} 


(5.1) 


(5.2) 


up  to  first  order  in  (x/R) .  R  is  the  radius  of  the  equilibrium 
Betatron  orbit.  Equations  (5.1)  are  non-linear.  To  use  perturb¬ 
ation  technique  let  us  first  write  them  in  the  following  way: 

x  +  ny^-fi^y 
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Writing 


where 


v  „  v(l)  .  v(2)  v  .  v(l)  .  (2) 

A  A  t  A  T  «  •  •  }  jr  jr  +  jr  t  •  •  # 


x(1)  +  nzy(1)  -  0 


J(2>  -  o  i«>  -  0 
*1 


with  the  given  initial  conditions,  and 


^)  +  nzy<2>«-nz  ^ 


.  (2) 


.(!)•(!) 


y(2)  -nzi 


(Did) 


(2)  .  n  /x-x\(2)  _  r  x(1)x 
UZ  l  R  /  UZ  R 


(5.3) 


(5.4) 


with  the  initial  conditions  x^2^  -  y^2^  •  0,  x^2^  -  y^2^  "  °* 
Higher  order  terms  will  be  neglected,  consistent  with  Eq.  (5.2). 

The  solution  to  Eq.  (5.3)  can  be  written  without  any  loss 

of  generality  as: 


x(1>  -  xrt  cos  fit 
o  z 

y<1)  '  ^sinnzt 


(5.5) 


where 


yG  “  flz  xo 


Inserting  this  into  Eq.  (5.4)  we  get  the  inhomogeneous  equa¬ 


tions 


X<2>  +  nz  y  (2^  -  -  nz  (i  +  cos  2  nzt) 


*(2)  -  nz  x  (2)  -  nz2  -fa  Bin2fizt 


(5.6) 


Thus,  in  addition  to  the  small  amplitude  oscillation,  there 


appears  a  drift  term  with  speed 

2 

Ymv 

_  X  C  _  c 

yD  2B  qR  “  ^  qR 


(5.9) 


where  u  is  the  electron  magnetic  moment. 

In  conclusion,  up  to  second  order  we  found  that  the  electron 
rotates  around  its  guiding  center  which  drifts  in  the  -y  direction 
with  constant  speed. 

Now  let  us  return  to  the  modified  Betatron  case.  For  simpli¬ 
city  assume  s  *  1/2.  The  equations  of  motion  up  to  second  order 
are: 


X  +  nz  y  +  UJ2X  -  -  nz  ^  y 
y  -  nz  i  +  U)2y  -  nz  |  i 


(5.10) 


where 


uu 


2 


(5.10') 


The  first  order  solution  is  given  by  a  combination  of  the  modes: 


U) 


± 


- 2 - 


(5.11) 


The  general  form  of  the  solution  is: 
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(5-12) 

where  the  constants  A,  B,  C,  D,  should  be  determined  by  the  four 
initial  conditions. 

It  will  be  shown  that  the  Fermi  drift  is  important  here  only 


when  uu  «  thus 


uj+  -  o 


(5.13) 


uu 


UT 

0. 


The  equations  for  the  second  order  corrections  are: 

x<2>  ♦  0i*<2>  ♦  *V2>  -  -  nz  .  «,(t) 


y<2>  -  nz*<2)  +  »2  y<2>  -  nz  |(1)i  (1)  .  ,  <t) 


(5.14) 


The  solution  to  these  equations  can  be  obtained  by  using  the 
method  called  variation  of  constants. 

Given  the  set  of  equations 


y  -  A(t)y  +  g(t) 


(5.15) 


where  A  is  anmxn matrix,  y,  y  and  g  are  n-vectors.  Assume  $  is 
the  fundamental  matrix  of  y  -  A(t)y.  Then  the  function 


■T^V.i’s*  ■*  ■“  ■*  rf 


A 
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u(t)  -  4>(t)  \  S>”  (s)  g(s)  ds 


(5.16) 


is  the  unique  solution  of  (5.15)  with  the  initial  conditions 
u(o)  *  0.  From  (5.12)  and  (5.16)  we  obtain  this  solution: 


U+-U)  L  o 


gx(S)  [sinuu+(t-s)  -  sinuu“(t-s)]  ds 


+  ^  gy(s)  [cos  id+  (t-s)  -  cos  ui“  (t-s)]  ds  L 

O  J 

-  — - —  K  g  (s)  [cos  uti  (t-s)  -  COS  UJ+  (t-s)]  ds 

«>-«»“  li  x 


(5.17) 


t  -j 

+  ^  gy(s)  [sinuj+(t-s)  -  sin  uu”  (t-s)]  ds  L 


The  other  terms  are  analogous  to  the  oscillating  and  the 
small  shifting  terms  in  (5.8).  Moreover,  it  is  only  the  constant 
part  of  gx(s)  which  contributes  to  the  Fermi  drift.  It  can  be 
shown  by  a  straightforward  calculation  that 


(5.19) 


7  ' 
yD 


w-o3  Lnz<x° +  y«J '  "  R*(x° +  yJ 


+  2a,2  (x  y  -  y  i  )1  SIBJlI. 
O  o  *o  o'J 


t  2  #  Q 

The  dominant  term  in  the  sum  is  fi  (x  +  y  ),  thus  we  get 

z  o  ** 


.(2) 


2m: 


sin  u d  t 


uu 


(5.20) 


The  Fermi  drift  speed  is  multiplied  by  cosuu”t,  and  the 
electron  will  drift  only  a  finite  distance,  Ay,  up  and  down. 


Uy| 


2RO  I  uu* 


i  lid  2 

55  hi 


(5.21) 


If  | Ay |  «  a  the  Fermi  drift  can  be  neglected. 

i2 


Close  to  injection  |ur|  -  CT  -  and  y^l,  thus  the  above 
criterion  is 


or  02  «  2nnaRr  (5.22) 

ox  c 

where  n  is  the  electron  density  and  r  is  the  electron  classical  radius 
11  -3  2 

For  n  - 10  cm  ,  0  -  0.5,  a -5  cm,  R  *  100  cm,  inequality  (5.22) 

is  satisfied. 

2  2 

When  y  Is  large  >  Cl  and  the  Fermi  drift  can  be  still 
negligible  due  to  the  focusing  effect  of  the  Betatron  field.  The 
condition  for  this  is  (vx/V||)2  «  a/R,  but  since  vx  -  v°/y  it  can 


be  written  as: 


< v2/v,|  )2  «  (a/R)  y2  (5.23) 

n  2  2 

Assuming  v  =“  c  and  (v°/c)  =  J  (5.23)  yields  y  »  10. 

2  2 

A  problem  exists  only  when  =“  ft  .  This  occurs  at  y  13 
(see  2.10).  In  order  to  prevent  the  beam  from  displacing  to  the 
wall  a  rapid  change  of  y  should  take  place.  However,  a  simple 
calculation  shows  that  y  should  be  changed  from  ~  12  to  ~  14  in 
about  30  nsec  which  is  not  practical.  It  turns  out,  however,  that 
the  surrounding  metalic  walls  help  in  preventing  a  catastrophe. 

In  the  preceding  discussion,  we  assumed  that  the  beam  is  located 
at  the  center  of  the  torus  and  neglected  interaction  with  the  walls. 
A  displacement:  of  the  beam  will  induce  non- negligible  electric  image 
forces.  The  equations  of  motion  for  the  beam  center  can  be  obtained 
by  averaging  (5-10)  over  all  the  electrons.  Let  us  denote  the 
coordinates  of  the  beam  center  by  X,  Y  and  the  single  particle 
coordinates  by  x,y: 


n_ 


x  +  nzY  + 


(5.24) 


0. 


Y  -  fi  X  +  u>*Y  -  <xx> 
z  b  R 


where  w2  -  *n2  -  (-f)  fi2  . 


(5.25 ' 


U) 


r.  is  the  beam  radius  and  ft. 
b  b 


* 


The  average  terms  can  be  calculated  immediately. 


fix  is  the  x  coordinate  of  a  single  electron  relative  to  the  beam 
/  -2, 

center.  (fix  )  *  0  because  the  beam  x-dimension  does  not  change 
as  it  drifts. 

/  2\ 

<  v  ) 

<xy>  =  XY  +  (fix-fiy)  =■  XY  +  - 


where  <6x*6y>  was  estimated  exactly  as  (5.20).  The  equations  of 
motion  of  the  beam  center  are: 


•*  2  < v, 2> 

X+V  +  V-  'IX 

Y  -  Q  X  +  u»?Y  -  XX 
z  b 


(5.26) 


The  constant  term  on  the  r.h.s.  result  in  a  small  shifting 

of  the  beam  in  the  -x  direction  (outside)  unless  uu.  *■0.  if 

,  D2. 

tab  -  0  a  drift  in  the  y  direction  with  a  speed  ~  will  start. 

z 

In  comparing  (5.26)  to  (5.10)  it  is  evident  that  the  beam 
and  the  single  particle  drifts  occur  at  different  yCcompare 
equations  3.10)  and  2.10).  Thus  when  the  single  particle  drift 
starts  the  beam  is  stable  and  it  will  only  oscillate  with  a  small 
amplitude.  In  order  to  estimate  this  single  particle  drift  effect 
assume  that  at  t  -  0  the  beam  is  located  at  X-Y-0,  and  it  starts 
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In  our  standard  example,  this  inequality  is  obviously  satis¬ 
fied. 

As  y  increases  further,  the  beam  will  become  unstable  for 
uub  -  0  andjwill  quickly  drift  to  the  wall  (se'e  Section  3  for  other 
Instabilities  which  occur  at  this  region) . 
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VI.  RESONANCE  INSTABILITIES 

In  the  previous  discussion  we  assumed  that  all  the  fields 
were  ideal,  i.e.,  they  did  not  depend  upon  the  coordinate  z. 

Errors  in  these  fields  give  rise  to  new  types  of  phenomena,  known 
as  orbital  resonances  or  parametric  instabilities.  These  in¬ 
stabilities  might  very  easily  destroy  the  beam  quality,  and  it  is 
our  purpose  here  to  give  quantitative  restrictions  on  the  fields  in  order 
to  assure  good  beam  quality.  Ve  shall  treat  the  equations  of 
motion  only  in  the  linear  approximation.  Higher  orders  are 
important  only  when  the  beam  quality  is  so  poor  that  they  are 
of  no  interest  anyway. 

In  the  previous  discussion,  the  fields  were  expanded  around 

the  equilibrium  radius,  R,  and  the  expansion  parameters  were 

assumed  constant.  Now  we  let  them  be  functions  of  z.  If  the 

electrons  move  with  constant  v  ,  these  parameters  are  also  periodic 

z 

functions  of  time,  with  a  period  T,  (t  ■  When  T  is  an 

z 

integral  (sometimes  half  integral)  multiple  of  the  transverse 
oscillation  time,  a  resonance  occurs.  We  shall  examine  some  of 
the  interesting  cases  and  find  the  restrictions  they  impose  on 
the  external  fields. 

We  shall  start  with  perturbing  the  Betatron  field  By.  This 
necessarily  perturbs  the  equilibrium  radius  R.  Electrons  will  be 
deflected  in  the  x  direction.  In  certain  conditions  they  will 
resonate  with  the  6  B„?(z)  field.  We  shall  derive  these  conditions 
and  determine  their  influence  on  the  operation  of  the  high  current 
Betatron.  6 By(z)  induces  additional  component  6Bz(y)  (v  x B  -  0). 
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This  will,  however,  add  only  second  order  terms  to  the  equations 
of  motion  and  thus  will  be  neglected. 

Assuming  6  B  “0,  the  Betatron  field  will  be 


By  “Bo  (1+x)  +  6  Vz> 


B 


x 


B, 


sx 

R 


(6.1) 


and  the  equations  of  motion. 

x  +  ix  +  nzy  +  [(i-s)  ny2  - 

y  +  |y*fizi  +  [sfly2'  fi2] 

In  discussing  (6.2)  let  us  start  with  the  conventional  Betatron, 
i.e.,  without  toroidal  field.  We  shall  see  later  that  the  high 
current  Betatron  can  be  analyzed  by  a  simple  generalization.  The 
equation  of  motion  for  the  x  component  is: 


21  2  6By(z) 
0  J  X  ”  fly  R  -Z— 


(6.2) 


y  -  0 


x  +  ix+  U)2x  -  no2  R  €(  t)  (6.3) 

^2  c  o  i  2  2 

where  fic  ~  Cly  ■  -g-  ^  ,  u)  *  Jfly  ”  and  it  changes  slowly  in 

2tt 

time.  €( t)  is  periodic  in  t  with  period  T  -  ?r-  ,  and  thus  can  be 


decomposed  in  a  Fourier  series: 
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r"? 

*'  t*v 

*.*•*  .* 
■*.  »v 


i 

& 

P>: 


■S'; 

i 

"■S- 


»-> 


6(t)  (oncosnfiot  +  gnsinnfiot 


T 

orn  -  i|  \  €(  t )  cos  n  fiQt 


dt 


(6.4) 


T 

^  ^  €(t)  sin  n  nQt  dt 

o 

If  uu  and  y  are  constant  the  solution  to  (6.3)  with  the  given 
initial  conditions  is: 

x  l  r  2 

x  -  xQcos  uut  +  —  sin  uut  +  —  \  fiQ  R  €(s)  sin  uj  (t-s)  ds  (6.5) 


/.***. 
**.  ^  • 
t#/£) 

r->> 


Substituting  (6.4)  for  €(s),  the  integral  can  be  calculated 
explicitly.  A  special  case  is  for  u>  -  nfi  .  Then  we  find  for  the 
integral  [_crnsin  nfiQt cos  nOQt  It  +  E  (non  secular  terms). 

Thus  the  amplitude  grows  linearly  with  time. 

A  conventional  Betatron  operates  below  the  n«l  resonance. 

In  the  modified  Betatron,  we  start  with  very  large  n  and  go  down 
to  the  lowest  resonances  as  y  increases.  It  is  thus  necessary  to 
estimate  the  amplitude  growth  as  resonances  are  crossed. 

For  large  n  resonances  are  crossed  so  quickly  that  no  growth 
(or  at  least  a  very  small  growth)  in  the  amplitude  is  expected. 

When  the  time  Atn,  to  go  from  uu  -  (n+i)fiD  to  w  -  (n-$)fiQ  is  com¬ 
parable  to  T,  the  time  to  go  around  the  torus,  a  real  growth  in 
the  amplitude  might  occur.  Atn  is  given  by: 

*  Ltn  "  n0  (6,6) 


duu 

Jt 


u)-nn 


Take,  for  simplicity,  only  the  cos  n  r^t  term  of  the  Fourier 
series  and  denote 


6u>  =  uu  -  n  0 


uS  *  (U)  +  n  ft0)/2 


then  the  resonant  term  will  be 


no  R<*n 
o  n 


sin  6  uj  t 


sin  w  t 


(6.7) 


Thus,  the  amplitude  for  6uu*o  goes  like  (6uu)“  .  This  is  the 

reason  that  real  growth  occurs  only  very  close  to  resonance,  and 

each  resonance  can  be  taken  separately.  Thus,  let  us  look  only  on 

the  n  resonance  as  it  is  crossed. 

If  A*  It  «  1,  the  asymptotic  expansion  WKB  can  be  employed. 
U) 

It  can  be  shown  then  that  the  solution  to  (6.3)  is: 


X  -  — —  cos  uu  t  +  sin  uu  t  +  R  ^  Jffy  sin  -  qp(s)]ds 


where 


qp(t)  -  J  u>(t')  dt' 


(6.8) 


In  the  high  current  Betatron  yuu  is  nearly  constant,  thus  the 
zero  order  amplitude  is  constant,  while  the  transverse  speed  de¬ 
creases  as  y"1.  The  integral  may,  however,  increase  the  amplitude 


due  to  the  resonance  terms,  which  can  be  shown  to  be 


2  It 

0.  Ran  c.  cos  [<r(s)-nQs] 


ds  sin  cp(  t) 


(6.9) 


The  last  integral  can  be  approximated  by: 


L*  vat/ 


UJ=nQ, 


2" 

s  ^ds 


2^2_  r 

j 


cos  t  dt 
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cos  t  dt 


(6.10) 


where  Nn  is  the  number  of  turns  the  electron  does  when  it  goes  from 

t  -  0  to  t~Atn.  We  assumed  that  Nn>l.  To  estimate  the  integral 

2 

(6.10)  we  note  that  the  function  cost  is  close  to  1  for  t<l  and 
is  highly  oscillating  for  t>l.  Thus  only  the  interval  o£ts*l  con¬ 
tributes  significantly  to  it.  Taking  33  -  nfi0,  the  integral  (6.10) 

is  estimated  to  be  —  and  substituting  back  in  (6.9)  we  get 

•Jrf  nn 

the  change  in  the  amplitude  as  we  cross  the  n  resonance: 


o_ 

£x  -  — s  yfr 

res  </*  n  n 


(6.11) 


To  be  specific  assume  that  the  error  is  constant  (AB)  in  an 


arc  L,  then: 


res  \  A 


(6.12) 


/  f 


33. 


To  minimize  AX„  _,  AB*L  has  to  be  smaller,  n  larger  and  the 

1*65 

time  to  cross  the  resonance  smaller. 

Now  let  us  return  to  the  modified  Betatron  case  in  details; 
The  equations  of  motion  are; 


V  •  •  2 

x  +  x  +  n  y  +  w  x 

T  ^ 


••  v  •  •  2 

y  +  ±  y  -  +  u)  y 

T 


R  e  (t) 


(6.13) 


The  change  of  the  solution  due  to  the  inhomogeneous  term  is 
given  by  the  integrals; 


Q?  R  J 


(t)  =  —2 - Ce(s)  [sin  uu+(t-s)  -  sin  uT(t-s)]  ds 

(u)+-uTW 


y(t> 


tip  K  r 
( uu+~uT )  ^ 


(6.14) 


€(s)  [cos  U) ” (t-s)  -  costu+(t-s)]  ds, 


where  uu^  are  given  by  Eq.  5.11.  Resonances  occur  for  u;*  -  mfi 

-2  c 
“  ^v  + 

However,  juu  |  “  <  thus  only  m  resonances  have  to  be 

wz 

considered  and  especially  the  low  n  resonances. 

Assuming  out  .  *■  2.1011  sec”1  (B^ -  104  Gauss),  and 
6  *1 

dy/dt  “■  10  sec  ,  we  obtain 


\  iSfi  -  h  Hfi  "  5l0‘6 « 1 


(6.15) 


and  WKB  can  be  employed.  The  results  of  the  one  dimensional  case 

+  2 

are  still  valid  except  for  an  additional  factor  ~(no/uu  )  .  Moreover, 
the  resonance  is  now  observed  also  in  the  y  direction.  In  our 


standard  case  (nQ/u)  )  «1,  almost  down  to  the  lowest  resonance, 


thus  reducing  the  growth  rates  predicted  earlier.  For  example, 

2 

if  n-2,  y  “*  340,  B„  =“  5  KGauss  and  (ft  /%  )  =*■  0.25. 


To  be  specific,  let  us  assume  that  ABQ  is  constant  in  an  arc 


L,  then: 


AX 


/ nx  .ab  n  .n  n2  /n“> 

(n)  *  (o\  L  (_o\ 

res  '  B  /  \  ,+/  n 


UU 


But 


Cl 


N. 


<_L  (*L\ 


U}+  T  ■ _  . _ 

n  1  Y-  Vdt/ 
n  n 


Thus  we  get  the  total  growth  in  the  amplitude: 


,AB 


i 


“res  '  hi2)  L  &)  2  1/0 


(6.1 


If  dy/dt  -  10®  sec”1  and  Clz°  -  2.1011  sec”1,  AXreg  *  15o(-g-2) 

o 

Most  of  this  comes  from  the  n  - 1  resonance  which  takes  a  very  long 

time  to  cross.  Thus  it  will  be  useful  either  to  cross  the  lowest 

resonances  very  fast  (by  turning  off  B_  in  t  «  10  jisec)  or  to  ex- 

z 

tract  the  beam  while  n>  1.  (e.g.,  by  using  the  drift  motion  of  the 

beam  -  see  Sections  ill  and  V). 

Before  discussing  the  general  case  let  us  work  out  another 

simple  case.  Assume  that  B_  is  perturbed  along  the  torus.  Accord- 

z 

ing  to  v  •  B  -  0,  ^  6BX  +  yp  6By  *  0  and  additional  first  order 
terms  should  be  included  in  the  equations  of  motion.  If 
yj  8BX  -  6By  and  if  6Bx(o,o)  »  6By(o,o)  -  0  the  equations  of 


motion  are: 


x  +  i  x  +  nzy  +  |  nzy  +  [|  n  2  -  n2]  x  -  o 


(6 


y  4  V  y  "  nzX  “  2 


1  a  „  n  -  2 

n_x  +  7T 


L2  “y 


n„  - Q  J  y-o 


Neglecting  the  terms  with  Y/Y  and  using  the  transformation 


x  +  iy 


exp  [i/2  J  n2(t')dt']  u ( t ) 


we  get 


u  +  u>  (z)  u  -  0 


where 


«  -  5  /O*  +  *0  *  -  40* 


When  uu  is  constant  the  general  solution  is: 


u  -  A  cos  uit  +  B  sin  uut  -  u  +  iu 

x  y 


and  u  ■  Reu,  u  ■  imu.  The  four  initial  conditions  are: 

*  y 


u  (o) 

X 


x(o) 


Uy(o) 


■  y(o) 


ux<°> 


a. 


n. 


(6 


(6, 


(6. 


(6. 


-  x(o)  +-g-  y(o) 


uy(o) 


-  y(o)  -  x(o) 


Unlike  Eqs.  (6.17),  (6.19)  for  u  and  u  are  not  coupled,  even 

x  y 


when  w  is  not  constant. 


For  the  stability  analysis  it  is  enough  to  examine 
|u(  -  |x  +  iy|.  The  amplitude  growth  rate  is  determined  by  the 
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two  real  equations 

..  2 

u.  +  tti  (z)  u.  -  0 
i  l 

where 

i  -  x,  y  . 

Hie  two  independent  solutions  to  each  of  these  equations  can 
be  written  according  to  Floquet's  theorem  as: 

u<t)  -  ei^it  g^(t)  (6.22) 

where  g^Ct+T)  -  g^(t)  and  the  constant  is  either  real  or  pure 
imaginary.  For  u>  constant  it  is  readily  verified  that  ■  ± 
However,  when  a  small  perturbation  is  added  to  uuo,  there  will 
be  regions  of  u>q  for  which  ^  will  become  pure  imaginary  and  one 
of  the  two  independent  solutions  will  be  unstable. 

To  show  this  explicitly  let  us  take  the  following  example, 
which  can  be  solved  analytically: 

Assume 

!uuo(1+Q)  0  <  t  <  a  -  T/M 

(6.23) 

<jjQ  a  <  t  <  T 

b  -  T  -  a,  and  uu(t  +  T)  ■  w(t).  It  can  be  shown  that: 

2 

cosuT  -  cos[u>0(T+  aQ)  j  -  2(^q)  sin  [u»Qa(  1  ♦  Q) j  sin  (uMj)  (6.24) 

As  Q  -•  0  the  r.h.s.  becomes  cos  u)QT.  Resonances  are  expected 
around  \uQT  -  nir .  if  Q  *  0  the  r.h.s.  might  be  greater  than  1 
(in  absolute  value)  in  some  intervals.  For  these  intervals  ti 
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will  be  imaginary  and  the  orbital  amplitude  will  grow  indefinitely 
in  time. 


In  order  to  find  the  locations  and  widths  of  these  intervals 
define 


X  «  ujq  [T  +  aQ] 


(6.25) 


The  r.h.s.  of  Eq.  (6.24)  can  be  written  as: 


(r.h.s. ) 


O  r 

cos  X  (1  +  G  sin  uj  b)  -  *  sin  2uu  b  sin  X 

o  £  o 


(6.26) 


where 


X _ 

H 1+0) 


If  |G|  «  1  the  r.h.s.  can  be  ±  1  only  when  cosX  =■  ±  1  or  Xs-  ntr, 


2fio 


(6.27) 


Expanding  the  r.h.s.  around  uuo  -  ur  '  we  get: 


2 

r.h.s.  -  (-)nCL+G  sin2tu(0n)b)  -  |sin2uu<n)b  (6x)  -|(6x)  ]  (6.28) 


The  maximum  of  this  expression  is  obtained  at 


6X  -  -  §sin  2  u,<n>  b 


and  is  1  +  Gsin2u/n^b  to  first  order  in  G.  The  width  of  the  n 


unstable  region  is 


fiUJ(n)  a.  sM  j8inuu<>n>b|  -  |Sin^n)b 


(6.29) 


V'  "V  V  -w  '  * 


For  small  n  resonances  (n«M) 


t  u,<n)  -  «  S  n 

o  M  o 


and  the  e-folding  time,  t^,  at  the  maximum  growth  rate  is 


(6.30) 


t  -  _L  .  M  .  T 
n  ttQ  n 


(6.31) 


In  the  high  current  Betatron  uaQ  changes  in  time,  due  to  the  change 

in  y  or  B_.  At  the  beginning  B„  is  held  constant  (~10kG)  and  y 

rises.  In  our  standard  example  ny  =**  670  and  unless  y  is  very  high 

n  »  1.  For  n  >  M  the  number  of  turns  that  an  electron 

4  O 

executes  while  it  crosses  the  n-th  resonance  is  £  10 

n 

Thus  until  n  is  very  low  this  number  is  much  smaller  than  1  and 

no  real  growth  is  expected.  As  n  becomes  of  the  order 

of  1  (n«M)  the  number  of  turns  the  electron  stays  in  the  n-th 

4  O  -2 

resonance  is  ~  3.10  — jj.  Assuming  Q  -  10  ,  M  -  10,  and  n  -  1  this 

number  is  ~  30.  The  amplitude  growth  will  be  less  than  10%.  Com¬ 
paring  to  Eq.  (6.16)  it  is  obvious  that  the  same  perturbation  in 
6By  is  much  more  destructive. 

As  we  pointed  out  earlier  the  missing  of  the  half  integral 

resonances  in  the  last  case  was  due  to  the  high  symmetry  of  the 

problem.  In  order  to  see  the  general  nature  of  the  motion  as  B_ 

z 

is  perturbed  assume  that  55B_/5x  -  0.  Thus  56  B  /5y  -  -  56  B  /5z 

x  y  z 

and  the  equations  of  motion  become 


x  +  y)  ■  0 

(6.32) 

y  -  fi  x  -  0 
z 


L-V-.’  v''  ■  V’1  v*  V' ’  V*'. ^  ’■  -.'Vi  k^L 


A •  V  -A'  .'-.  A. 
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Integrating  the  first  equation  and  substituting  in  the  second  we 
get 


x  +  Q  y  =  +  (n  y)^°^  =  A 

Z  Z 

y  +  nz2y  =  a  nz  -  A(nz  +  6nz> 


(6.33) 


A 

where  0  is  the  average  along  the  torus  and  A  is  a  constant  which 

Z 

depends  upon  the  initial  conditions.  The  second  equation  is  a 
nonhomogeneous  Hill’s  equation  and  it  has  integral  and  half  inte¬ 
gral  resonances.  The  symmetry  in  this  case  is  obviously  lower  than 
in  the  previous  case. 


Let  us  now  consider  the  most  general  perturbations  in  the 
fields  for  a  high  current  Betatron.  These  fields  may  be  expanded 
around  the  average  equilibrium  radius  R  -  vz/fiQ.  Using  Maxwell's 
equations  v  •  B  «  0,  7  x  B  ■  0  and  keeping  terms  only  up  to  first 
order  in  the  equations  of  motion  we  get: 


x  +  i  *  +  o  y  +  <i  2  -  n2)  x 
Y  ^  ^  y 


n  6B  9  /6B  \ 

ny  R-g-  +  ny  ^  +  68  : 


y  +  ( 1  -  or) 


y]  =  f(t) 


[•■51 
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••  y  •  •  X  ^ 

y  +  -h  y  -  ftzx  4  n  -  c  )  y 


6BZ  •  SBZ  1 

+  nz  LT"  x  +  a  —  XJ  ‘  E(t> 
z  z 


(6.34) 


where  we  assumed  s  =  1/2  and: 


Bz  -  Bz  +  6Bz(z) 


By  ■  B0(l+f)  +  tBo(z)  +  ^+  M^)bo  f  +  («.(«, -X) 


3  6B„ 


oz 


(6.35) 


SB*<z> +  (^f +  ¥)  Bo  ¥  -  “(z) 


3  6B2 
dz 


x 


a(z)  is  defined  by: 


36B? 

“ST 


-  o(z) 


S6Bz 

“Si”  • 


We  wrote  the  equations  of  motion  so  that  all  the  perturba¬ 
tion  terms  appear  on  the  r.h.s.  These  terms  are,  however,  not  of 
the  same  order.  The  dominant  terms  are  the  first,  thus  one  sus¬ 
pects  them  to  be  the  most  dangerous  for  the  beam  quality.  If  6BQ/B 
and  6B  /B  are  given,  their  contribution  can  be  calculated  accord- 

X  X 

ing  to  Eq.  (5.17) . 


The  first  terms  result  in  resonances  for  u>  -  nn  (n  >C1  ). 

o  z  y 


These  resonances  were  considered  in  detail  (see  Eqs.  (6.1)  -  (6.16) 


The  only  difference  is  that  here  B  (z)  is  perturbed  too. 

A 


The  growth  rate  at  the  resonance  is  determined  mainly  by  the 


larger  of  6B  and  6Br.  Thus  to  minimize  the  effect  both  6B  and 

X  j  X 


6By  should  be  decreased. 


The  second  largest  terms  on  the  r.h.s.  of  Eqs.  (6.34)  for 


n  >  o  are  the  n  terms.  We  have  seen  earlier  (see  Eqs.  (6.17)  - 
Z  y  z 


(6.31))  that  for  a  *  1/2  these  terms  give  rise  to  parametric  in- 

j — 5 - 2 - 5" 

stabilities  with  yn  +  2fi  -40  •  n0  .  For  a  *  1/2  new  re¬ 


sonances  appear  (|u>  |  -  1/2  n0Q,  Gz  -  nCJQ) .  These  new  modes  disappea 


as  a  -  1/2.  As  III"  -  0  and  u>+  -  0_  (only  B_  is  left)  we  recover  the 


-  -S  fi  resonances  (see  Eqs.  (6.32,(6.33)). 
z  z  o 


The  other  terms (proportional  to  6BQ/2Bo  +  6s)  give  rise  to 


resonances  at  uu  -  nn 


All  the  instabilities  we  have  just  considered  start  growing 
linearly.  Comparing  the  terms  of  the  r.h.s.  of  Eqs.  (6.34)  it 
is  obvious  that  the  first  terms  are  greater  by  ~  R/rb  than  the 
second.  This  makes  the  first  terms  much  more  destructive.  The 


remaining  0  terms  were  investigated  carefully  and  it  was  found 
z 


that  they  are  at  least  one  order  of  magnitude  less  dangerous  than 


the  first  terms. 


DISCUSSION 


In  this  paper  we  investigated  some  aspects  of  the  stability 
of  a  High  Current  Betatron.  Specifically,  ybg  discussed  (1)  space 
charge  limits  of  both  single  particle  and  the  electron  ring,  (2) 
the  Fermi  drift,  (3)  resonance  instabilities  and  (4)  beam  emitt- 
ance.  '  _ _ 

Two  important  questions  have  not  been  addressed,  namely; 

(i)  stability  during  initial  acceleration  (injection,  trapping 
and  initiating  acceleration  will  be  discussed  separately),  (ii) 
collecrive  modes,  such  as  negative  mass  and  resistive  walls. 

These  modes  have  been  discussed  extensively  by  many  authors,  yet 
there  is  no  conclusive  evidence  that  they  might  cause  a  destruc¬ 
tive  damage  to  the  High  Current  Betatron. 

Let  us  summarize  the  main  results: 

(1)  The  space  charge  limit  during  injection  is  increased 

2  2 

compared  to  a  conventional  Betatron  by  (B_  /2b„  ) . 

z  y 

(2)  Single  particle  orbits  become  unstable  (for  s  * 1/2)  at 

3  2  2 

Y  «  4NR  rc/rb  •  This  instability,  however,  is  self-stabilized  re¬ 
sulting  in  only  a  small  increase  of  the  beam  radius. 

(3)  The  stability  of  the  beam  is  influenced  by  the  surround¬ 
ing  walls.  Ideal  conducting  walls  enhance  the  stability  of  the 
beam  relative  to  the  single  particle.  In  the  opposite  case,  when 
the  image  currents  decay  instantaneously  the  beam  may  become  un¬ 
stable  due  to  space  charge  limit  as  y  increases,  even  for  the 

s  »  1/2  case.  This  occurs  only  if  the  line  density  is  high  enough. 

2 

(4)  For  s  *  1/2  the  ring  becomes  unstable  at  y  ^  4N(R/a)  r 


(see  Eq.  (3.10)).  In  this  region  the  total  focusing  force  on 
the  beam  becomes  zero  and  various  phenomena  occur  which  make  it 
practically  impossible  to  hold  the  beam  in  the  torus. 

(5)  The  Fermi  drift  is  negligible  except  in  the  region 
where  the  total  focusing  forces  on  the  beam  become  zero  (this  is 
true  even  for  s  -  1/2).  Then,  the  beam  drifts  quickly  to  the 
wall.  This  drift  may  be  utilized  to  extract  the  beam. 

(6)  The  beam  emittance  changes  only  slightly  during  accel¬ 
eration  (neglecting  instabilities).  It  is,  thus,  determined  by 
the  injection  method  and  by  the  instabilities  enhancement. 

(7)  Errors  in  the  external  fields  may  cause  driven  as  well 
as  parametric  resonances.  The  first  are  caused  by  errors  in 

B  (z)  or  B  (z)  which  are  not  zero  at  the  equilibrium  orbit  (mono- 

y  * 

pole-like) .  These  are  the  most  dangerous  resonances  and  put 

severe  restrictions  on  those  fields  unless  the  beam  is  extracted 

before  they  develop  (n>l).  Parametric  resonances  (dipole-like), 

caused  by  errors  in  B  (z) ,  B  (z)  or  s(z)  are  much  less  destructive. 

z  o 
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TABLE  I:  Standard  set  of  parameters 
for  the  high  current  Betatron  (used 
throughout  this  paper). 


Major  radius 

(R) 

100 

cm 

Minor  radius 

(a) 

5 

cm 

Bz 

10 

KGauss 

Line  density 

(N) 

2.1011 

-1 

cm 

Beam  radius 

<V 

1 

cm 

2.1011 

sec  ^ 

•l 

1.4.1010 

sec-1 

bd/sy 

- 

16 

Gauss 

-1 

dY/dt 

1 

Msec 

Fig.  1. 


Figure  Captions 


Schematic  design  of  High  Current  Betatron  with 
inductive  charging. 

Fig.  2.  Coordinates  for  a  toroidal  electron  beam. 

Fig.  3.  Graphical  representation  of  the  stability  con¬ 
ditions  for  the  modified  Betatron  (see  Eq.  2.8)) 
The  transition  takes  place  when  Y  is  increased. 


